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They are the basic values for the computation of Mie scattering coefficients. The computations have 


been carried out with the IBM 701 electronic data processing machine. 
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i, The Probiem 


The scattering of light in the atmosphere has become of increasing interest and importance to 
research in certain fields of physics, meteorology and colloid chemistry. A sound knowledge of the 
scattering characteristics of particles of size comparable to the wavelength of the incident light is 


necessary in studying the effects of light scattering. 


This problem has been solved theoretically over forty-seven years ago by G. Mie!. His 
formulas are exact solutions of Maxwell’s equations for spherical particles of any size and material, 
lie integrated Maxwell’s equations for the electric and magnetic field inside and outside a particle and 
culcalared the amplitudes of the waves of light radiated by the particle in various directions. The exact 
sohitions appear in the form of series which converge more slowly with increasing radius of the particle 


if the waveiength of incident light remains fixed. 


Various investigatorshave since used his formulas to carry out numerical computations. Unfortu- 
nately, chis becomes an extremely time-consuming task for large radii due to the slow convergence of 
the series, Some of the earlier computations are unreliable for large radii because not enough terms 
were carried. As a rough approximation it may be said that the number of terms in the series which are 
required before truncation, is equal to the value of 4, the size parameter. Hence, the available data 


are limited. 


The only reliable results for large size parameters, that is, particles of large radius, are pub- 
lished by Gumprecht and Sliepcevich?.. However, they used very large increments of the size parameter. 
Yne existing tabies do not meet the need we have for such tables; therefore, we planned and carried out 


aew computations, 


The diffraction theory undertakes to solve Maxwell’s equation for a plane wave coming trom 
iafinicy hitting any surface of discontinuity of the optical parameter (index of refraction, n). An exact 
sohiticey can be found only for simple geometric surfaces of discontinuity, in such a case the equation 


reduces to a boundary value problem for ordinary differential equations. 


Mie! solved the case for an incident plane wave falling upon a sphere of radius r. The sur- 
rounting medium is transparent, homogeneous and isotropic. The incident wave induces forced oscil- 
laiioas of free and bound charges which are synchronous with the applied field, inducing a secondary 
clecuic and magnetic field each of which consists of two parts, namely one inside the sphere and an- 
other om. curside the sphere. We ask for the intensity of the scattered light as a function of the scat- 
teeing argyle . ‘Thus, the intensity scattered, J, is a function of four parameters, r, A, 9, and n, namely. 


}=Jof {er A,0, Of. 21) 
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~ The solution has the form of infinite series, the coefficients of which are in turn determined 

by infinite series of Bessel functions and Legendre polynomials. 

“ We define the following terms: . 

ns n = index of refraction (of the aerosol particle against the dispersion medium) 

ne r = radius of spherical particle (aerosol particle) 
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h= wavelength of incident light in the surrounding medium 


ae 2h '/\ = size parameter 


s 


és 4 B= a G& 
“z = scattering angle (angle becween the incident and the scattered beam) 
x = COS ) 

A Plane of observation = plane containing the direction of the incident wave and the direction 
x? of the scattered wave. 

“ , 
ss if the incident cadiation is unpolarized (natural light) and its intensity is J, then the intensity of the 
ey radiation scattered in the direction 9 is given by Mie as 

"5 1 , : 
te j{o, 6} = Jo a2 4) {a, a, Of} + 19 {a, n, 8} ° (2.2) 
me ‘The quanticies +; sad ty ave called intensity functions, They have the following meaning: The intensity 
aot 


functions i, and i, are proportional to the magnitude of the two incoherent, plane polarized components 
scattered by one particle; i, is proportional to the intensity whose electric oscillations are perpendicular 
to the plane of observation, i.e, polarized light in the horizontal plane; i, is proportional to the intensity 


F 2 
ry 


oe 
—o whose electric osciliations are in the plane of observation, i.e. polarized in the vertical plane. 
ae This resule already shows that the scattered radiation is always polarized if the incident radiation 
i is unpolarized. The intensity functions* derived by Mie are 
oo 

; 2m+1 2 
me i, = 4 a b. T 
ou 1 re m(m+1) m “np * °m ‘m 
ol 
ee (2.3) 
oo 
tae 
" as 2 
me "275 > 2m +1 am To t bm “n 
wey m (m+) P 
on 1 
we 
hee 
ws The vertical fines indicate that the absolute value of the complex argument isto be used. It is seen that 
ey the scauered wave is composed of vatious partial waves and the amplitudes of these partial waves depend | 
sy ona, and b,,. Phe cerms in these series are defined as follows. The emplitude functions a,, and b,, are | 
ian complex numbers: ! 
aN ' 
me ae 
eS * The noration used here differs somewhat from the notation used by Mie! and those by Lowan. It also differs 
ws slightly from she first nate oa those studies written by one of us (B. Goldberg‘). The formulae are written in 


such a way 4s fo make wnmachine computations simpler. 
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1/2 : 
S, {z} = = (er 4 1/2) 
1/2 
Ca {z} = CD™ — - ee 


aa {z} = AG (2.5) 
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" 
~ Yhe argument z takes on the values Qor B in Eq. (2.4), and 
: | mti/2 = Bessel function of positive half order, 
a “ 
ee | “m-1/2 = Bessel function of negative half order. 
) ; 
ye The amplitude functions a,, and b,, can be interpreted physically as the m’ch electric parcul 
om wave and the m’th magnetic parcial wave, respectively, The solution (2.4) vives che scactcred field as 
va 4 spherical wave composed of two groups of partial waves. The first group, the electric partial wave 
" 4,, is that pace of the solucion for which the radial component of rhe wwagneric vector in the incident 
Be wave is zero. Similarly, the second group, the magnetic partia) wave by, refers to zero radia} compo- 
ce nent of the elecuic vector. 
a ‘the values of a, and b,, dependon the material (index of refraction n) and the size of the 
hoe particle (size parameter (). The partial waves can be considered as coming from an electric or 
= magnetic muitipole field, The first partial wave comes from a dipole, tie second from a quadrupole 
wwe and so oa. 
*r . 
ge The results show that the amplitude function Ang Of the (midith << strie partiat wave is ul 
ae the same oder of magnitude as b,, of the m’th magnetic partial wave as alieady mentioned by Mic 
sa 2 and can be checked by any table listing the values for ay, and h) Cfor exeapte, Lowanj.? An ex- 
ample may show this. For @ ~ 0.1 only 1 value of a, and by bas the go uteot order of niannituds 
“7 (107! and 10°?), which means for (= 1 the first partial wave dominates, aswdtfor @- 20 about 22 
t 
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partial waves have amplitudes of the same order of magnitude. The number of the partial waves - neces- 


sary for numerical compucations + increases with Q. In general, the number of terms in the series in- 
creases with 4, and above Q= 10 is given by 

My=7 + 1.2 a, (2.6) 
for our values computed for &=0 to 30 and the index of refraction in the range n = 1.33 to 1.50. For 


< 1G the number of terms is 


a=1 2 3 4 5 6 7 8 9 10 
m 5 7 9 10 12 13 14 16 17° 19. 


the coefficients 7, and T, are essentially angular functions containing the scattering angle @, buc 


are completely independent of a and a: 


dP d2 p 
m m 
T = T -—— 2.7 
- d x m d x2 ( ) 
= 2 
Tn =x Na - (1 — x*) Ty 


#2 2 cos? e-2 and 1 28 cos? 0 ~— 33 cos 8. 


we shail discuss these functions further in part 7 of the publication where we compute the angular Mie 
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Z tarioes of Aroplitade Functions a, and b,, 
Phe compurstioms for a,, and bas can be seen from equations (2.4) involve the calculations of 
OES Fes t ae { zy 4 Se { z} and C’, { z } for the values of m required for convergence 
OF Jroosuriva, in proyramming the computations, these functions were computed using recursion relation- 
SDays crtsanig tan those satisfied by the Bessel functions: 
‘ i) | sing 3; C = 
et : + Co Yzt =cosz , 
‘ {2 - Stn Zz . = cos Z . 
ane cy a ~ cos z ; Cy {z} as + sinz , 
( > Paria { } 2m-1 
peepee S = . “= - 
m bet Fa m= Sm-2 + Cm Uz Zz Cel Cm-2 (3.1) 
s! 4 ) 7 m % 
ail wee ee gs Sm + Smet ’ 
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Bi Then we write 
tee 
: += sy {Bh sy fa} - 0S (ah oe {8h 
ne 
“ es Ss). {B} Ca {a} - aS, iB}, {a} : (33 
" -3) 
cE u=nsty {B}s, {a} - sty fe} sn d8}. 
oN Neva 10 Sn {8} Cm {a} ae oa fa} Sm {3} } 
oS The amplitude functions Eq. (2.4) can now be written as 
A 3 _ m+1/2 s _ m+1/2 Ven tay Se Pots 
y fa 18) CD See ed Tear FT a] 
» or 
3 m 1 t/s 
> ay {af = CD eee! 1 eee 
= m [ 7 a ee | a2) 
a ( 1+ (/sye 1 + (/s)? , 
A similar derivation yields b,.. Hence the amplitude functions are givea by 
; _ -D™ (t/s) 50 Is cc ] 
am fan} = 4 r/se Fh Le (se? aa 
Ak 
(-1)mt1(v/u)  , (-a mt \ 
b a, = ———_— | a ue 
a a} 1 + (v/u)2 Be 1 + (v/u)2 
or in the standard terminology : 
( 
a, {%nf = Ra) + Ta), | 
‘ 7 CSE? 
bm {anf = R (by) + F (by) 


In computing the terms of the series, it is found that as Qincreases the number of terms required tue 
convecpeace increases. The terms oscillate about +] up to some m and then bein to decrease very rapidly, 
In the recursion relationships for S,, and C,, for small z, accuracy is lost fairly quickly as m increases, bus 
since the number of terms required is small, this is of no consequence. For large z, where the number of 
terms required is increaSed, there is again no difficulty because of the division by z. lt was alse found that 
u and v, not much accuracy is lost duc to the subtractions until the terms of the series be- 


ia computing s, 4, 


come sufficiently small. When this is the case the number of digits required in that teem ts not as great. 
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The derivatives were obtained using the previous functions rather than the previous derivatives so that the 
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error due to differentiation was made as small as possible. Our definitinn o- a, and ba (2.42 dirtors rions 


ws 


Taya 


that used by Lawan aad by Gumpreche and Sliepcevich by a factor 


2m + 1 


Yhis faccot aupears in our definirion of iy and iy (2.3), since itis more advantageous from the com- 


putatiosal scindpoint, Our ty and iy values agree numerically with those in the other tables. 


These computations were made using the IBM 701 Electronic Data Processing Machine and were 
pre grasa] asin, UM's floating poiat system, Speedcoding I. The sine and cosine subroutines used in 
the compa i nafy cos OS. and C, give values which are good in the range considered, to eight 
Sev val adares. Phsiatecance placed on the truncation of the terms was 10°8 for both Real (a,,) and 


‘+ a, and b, are believed to be valid to 6 significant figures. The results have 
i. Jans? existing tables, like Lowan’s® and Gumprecht and Sliepcevich’s?. The results 


tooo ota tp other tables. 
neve toatisi machine ecrors, every calculation was duplicated and a check sum carried with 
ciick sums were matched by machine and the calculation rejected on non-match. 
‘oy se oes ut transeription were eliminated by machine-printing these tables by the IBM 701 
‘.. resulting master sheets. The results were printed outto 10 digits, but the last 


Pee wr gut out hefore the offset printing. 


ranged for increasing values of &. Under each & value all terms of the amplitude 


ui all four functions R (a,), I (a,), R (b,,) and I (b,,) are smaller than 1 x 


re “os 9 teras Fulfill this condition. But for & = 30, 43 terms are necessary. 
Sos anne list. the number of the term. The second column gives R(a,,): the third K(a,,), 
ne . omtsie iant (bp) ). The numbers are printed to read 0.xxx, i.e. the decimal point 


‘a, fies fipure. [te is followed by the exponent to the base 10. They are printed 
Seuelelie = yo 
ee, 


sige 7 should read: 
Rfa,) = + 0.1914571 x 1079 
K(a,,) = - 0.3665583 x 107 
R(b,,) = - 0.2663852 x 10° 
(b,.) = + 0.7096110 x 10°)3 
Pe ca ts yrreGi en dipis. The last three were covered before offset printing. Hence, 
be fetes ai st tn iydats, the last one is not rounded off. Six digits are believed to be valid and the 


senrabh ara sols bte so that the reader can round off, if he wants six digits for further com~ 
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